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Abstract 

We consider the effect of weak uncorrelated quenched disorder 
(point defects) on a strongly fluctuating flux-line liquid. We use a 
hydrodynamic model which is based on mapping the flux-line system 
onto a quantum liquid of relativistic charged bosons in 2+1 dimensions 
[P. Benetatos and M. C. Marchetti, Phys. Rev. B 64, 054518, (2001)]. 
In this model, flux lines are allowed to be arbitrarily curved and can 
even form closed loops. Point defects can be scalar or polar. In the 
latter case, the direction of their dipole moments can be random or 
correlated. Within the Gaussian approximation of our hydrodynamic 
model, we calculate disorder-induced corrections to the correlation 
functions of the flux-line fields and the elastic moduli of the flux-line 
liquid. We find that scalar disorder enhances loop nucleation, and 
polar (magnetic) defects decrease the tilt modulus. 
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1 Introduction 



The physics of vortex-line arrays in type-II superconductors has been the 
subject of intense research activity in the past fifteen odd years [1, 2, 3, 4, 
5] . Potential technological applications of high-temperature superconductors 
rely on the pinning of the vortices in order to eliminate dissipative losses from 
their motion. In addition, the competition between thermal fluctuations, 
intervortex interactions, and various types of disorder gives rise to a variety 
of phases, thus making vortex matter an excellent laboratory for the study 
of fundamental statistical mechanics. 

The effect of disorder on the Abrikosov lattice phase of vortex arrays has 
been the subject of extensive theoretical and experimental studies [6, 7, 8, 9, 
10, 11]. Much less understood is the role of disorder in the vortex liquid phase, 
which is known to occupy a broad region of the magnetic-field-temperature 
phase diagram in high-temperature superconductors. A powerful method for 
describing both clean and disordered vortex liquids is the mapping of the 
statistical mechanics of directed line liquids in three dimensions onto that of 
two-dimensional nonrelativistic quantum bosons [12, 13, 14]. This mapping 
assumes small fluctuations of the flux lines away from the direction z of the 
external applied field and intervortex interactions only among line segments 
at the same "height" , z. The vortex lines map onto to the world lines of the 
quantum particles and this mapping is particularly useful for describing the 
effect of columnar pins caused by heavy-ion irradiation, which correspond 
to quenched impurities for the bosons [15]. In contrast, point defects in 
the superconductor correspond to time- dependent disorder for the bosons 
and are more difficult to treat in the quantum system. It has been shown, 
however, that this type of disorder yields "Lorentzian-squared" corrections 
to the thermal static structure function of the vortex liquid [16, 17] and does 
not renormalize the tilt modulus in the thermodynamic limit [17]. 

A more convenient framework for describing point disorder is the hydro- 
dynamic description of directed flux-line liquids developed some time ago 
by M. C. Marchetti and D. R. Nelson [18, 19]. This earlier implementa- 
tion of vortex liquid hydrodynamics takes into account the repulsion among 
flux-line segments at different heights, z, which renders the intervortex in- 
teraction nonlocal, but is restricted to directed lines, excluding overhangs or 
spontaneous vortex loops. 

It was argued recently that vortex loops may be important in the vor- 
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tex liquid phase, especially at low fields, where diverging loop fluctuations 
may destroy the coherence of the vortex state via a so-called "loop blow-out" 
transition [20, 21, 22]. Although earlier numerical evidence for such a tran- 
sition in the uniformly frustrated XY model [22] was subsequently shown 
to be associated with boundary effects [23], these simulations do indicate 
that spontaneous vortex loop formation is important in vortex liquids. Also, 
recent experiments in YBCO show evidence of a liquid-to-liquid transition 
which appears to agree with the "loop blow-out" scenario [24]. Disorder, 
which is ubiquitous in real high-temperature materials, was not included in 
the above mentioned simulations. This renders a detailed comparison with 
the experiment problematic. 

In this paper we present a hydrodynamic model of strongly fluctuating 
flux liquids that allows us to describe in a unified manner both field-induced 
lines and spontaneous loops. The key ingredient is a new scalar field that 
describes the density of flux-line length regardless of its direction or orienta- 
tion and therefore incorporates the contribution from neighboring antiparallel 
segments associated with closed loos. We use this framework to study the 
effect of point disorder on the flux-line liquid by evaluating perturbatively 
disorder-induced corrections to the static correlation functions and to the 
elastic constants. Both "scalar" point disorder (that couples to flux lines 
irrespective of their orientation) and "vector" point disorder (arising from 
impurities carrying magnetic moments) can naturally be introduced in the 
model and play qualitatively different roles in renormalizing the elastic prop- 
erties of the liquid. In particular, "scalar" disorder softens the single-line 
energy, enhancing the formation of vortex loops. 

The statistical mechanics of strongly fluctuating vortex liquids in zero 
external magnetic field has been described theoretically by M. Kiometzis et 
al. [25] and by Z. Tesanovic [20, 21] by mapping the vortex line system onto 
a collection of relativistic two-dimensional quantum bosons. In this map- 
ping spontaneous vortex loop formation corresponds to particle-antiparticle 
creation and annihilation. Building on this work and its generalization to 
vortex liquids in a field, we recently proposed a general mapping of flux-line 
liquids of directed field-induced lines and spontaneous flux loops onto a two- 
dimensional system of relativistic charged quantum bosons, where the exter- 
nal applied field for the bosons corresponds to the bosonic chemical potential 
[26]. The main difference between our model and earlier ones [20, 21, 22] is 
that, while these authors consider fluctuating vortices in a frozen magnetic 
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field (corresponding to the limit of infinite penetration depth), we describe 
magnetic flux lines fluctuating in concert with their vortex cores by keeping 
the London penetration depth finite. By coarse-graining the boson model, 
we obtained a generalized hydrodynamic description of strongly fluctuating 
liquids of flux lines and loops. 

After describing the model in Section 2, we use a Gaussian approximation 
of our generalized hydrodynamic free energy to evaluate the renormalization 
of the liquid properties from three different types of uncorrelated disorder in 
Sections 3, 4 and 5. Finally, we summarize in Section 6. 

2 The model 

In the hydrodynamic description of a directed flux-line liquid proposed by 
Marchetti and Nelson [18, 19], the large scale properties of the system are 
described in terms of the three components of a coarse-grained vector field, 
t(r) = {t±,t z }. The component t z is the areal density of flux lines along the 
direction z of the external field piercing the xy plane and directed along z. It 
is proportional to the z component of the local magnetic field in the material. 
The two-dimensional vector field tj_ describes the local tilt away from the z 
direction and is proportional to the local magnetic field in the xy plane. The 
correspondence between this hydrodynamic model and the familiar nonrel- 
ativistic boson model is well understood. It is important to note that the 
energy of a collection of flux lines in an applied field has full three-dimensional 
rotational symmetry. This symmetry is reduced in the nonrelativistic boson 
model which is invariant under Galilean boosts, corresponding to a uniform 
tilt of the flux lines away from the direction of the applied field. In the con- 
ventional hydrodynamic free energy, the single-line part is invariant under 
uniform tilt (Galilean symmetry), while the intervortex interaction has full 
rotational invariance in three dimensions. Full rotational invariance of the 
flux-line system corresponds to Lorentz invariance of the bosons, indicating 
that the quantum particles should be treated relativistically. 

In Ref. [26], we showed that a three-dimensional liquid of directed field- 
induced lines and spontaneous flux loops can be mapped onto a (2 + 1)- 
dimensional liquid of relativistic charged bosons interacting via a screened 
electromagnetic field. We exploited this mapping to obtain a general hy- 
drodynamic description of strongly fluctuating vortex liquids that allows for 
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overhangs, loop formation, and nonlocal (in z) interaction among vortex 
segments at different "heights". The hydrodynamic description of strongly 
fluctuating vortex liquids is in terms of four coarse-grained fields: the three 
components of the vector field, t(r) = {tj_, t z }, which measures the local mag- 
netic field in the material, and an additional scalar field, p(r). The new scalar 
field p(r) appears because of the possibility of spontaneous loop formation. 
It measures the total contour length of flux lines enclosed in a unit volume 
centered at point r, regardless of their direction or orientation. Upon coarse- 
graining, configurations of antiparallel neighboring flux lines, as arising from 
the spontaneous formation of a vortex loop, do not contribute to t(r), but do 
contribute to p(r). In the relativistic boson mapping, t(r) = {tj_(r), t z (r)} 
corresponds to the conserved current- c/iarge-density, which is a (2 + l)-vector 
with the time-like component t z the conventional charge density of bosons 
and tj_ the charge-current density. In contrast to the charge density, t z , the 
number density of bosons is not conserved due to creation-annihilation events 
(corresponding to vortex loop formation in the flux-line system). The scalar 
field p(r) corresponds to the square of the amplitude of the complex boson 

field, $(r) = \J p(r) exp[i0(r)]. The correlator ($*(r')$(r)) is proportional to 
the probability of finding a flux line segment connecting the two points, r 
and r' [27], irrespective of the segment's orientation. 

The hydrodynamic free energy functional of a strongly fluctuating flux- 
line liquid with disorder can be written as T = J-q + J-'d, where JF is the free 
energy of the pure system and Tr> incorporates the coupling to point defects. 
The pure contribution for a line liquid in an external field, H = H z, is [26] 

+^ i P'(?>|t(q)l 2 , (1) 

where e\ = eoln/t is the bare flux-line tension, eo = <Pq/ '(47rA) 2 , A is the 
effective London penetration depth, and k = A/£ is the Ginzburg-Landau 
parameter, with £ the coherence length. The last term (written in Fourier 
space) describes the intervortex interaction, with V(q) = 47re A 2 /(l + g 2 A 2 ). 
The average equilibrium values of the hydrodynamic fields are: tj_o = 0, 
t z0 = ±p , and p = (\H \ - H cl )/<f> for \H \ > H cl , while p = for 
\Hq I < H cl . Here H cl = 47rei/0 o is the lower critical field and p = 1/oq 
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represents the average areal density of flux lines in the xy plane, with ao 
the intervortex spacing. Averages over the conformations of flux lines and 
flux loops are evaluated with the Boltzmann weight ~ exp[— jF//c B T] and the 
constraint 

V-t = 0, (2) 

which prohibits flux lines from starting or ending within the sample. In the 
relativistic boson mapping it corresponds to charge conservation. In this 
paper we present the hydrodynamics of vortex liquids in isotropic super- 
conductors, but the generalization to the anisotropic case is straightforward 
[26]. 

In contrast with the hydrodynamics of directed flux-line liquids, the ef- 
fective free energy of Eq. (1) respects the rotational invariance of a pure 
flux-line system in an external field. This rotational invariance corresponds 
to the Lorentz invariance of the relativistic boson liquid. The non-Gaussian 
hydrodynamic free energy of directed vortex liquids [28] is obtained from Eq. 
(1) simply by setting p(r) = t z (r). This suggests that the difference 

tf(r) = t z (r) - p(r) (3) 

may be interpreted as a density of lines associated with spontaneous vortex 
loop fluctuations. 



3 Point defects with magnetic moments par- 
allel to the external field 

Uncorrelated quenched disorder from point material defects is usually mod- 
eled by a random field which couples to the fluctuations in t z (r) [16]. The 
disorder contribution to the free energy is written as 

Td = J d 2 r ± dzV D (r)5t z (r) . (4) 

The quenched random potential is taken to be Gaussian, statistically homo- 
geneous and isotropic, with Vb(r) = and Vd^)!^^) = A5(r — r'), where 
the overbar denotes the quenched average over disorder. A random field that 
couples to t z corresponds to quenched magnetic dipoles of random strength 
and position, with magnetic moments aligned with the z direction. It favors 
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local alignment of the flux lines with the direction of the external field (when 
Vo(r) < 0). It corresponds to a random chemical potential for the nonrela- 
tivistic bosons [17], and to a random scalar potential coupling to the charge 
density for relativistic charged bosons [26]. 

In order to calculate the disorder-induced contributions to the various 
correlation functions characterizing the vortex liquid, we expand the flux- 
liquid free energy of Eq. (1) to quadratic order in the fluctuations of the 
fields from their equilibrium values, 5t(r) = t(r) — t and 5p(r) = p(r) — p . 
The resulting Gaussian approximation is given by 

F G [5t,5 P ] * ^ + ^E([^ + %v ! ^]lMq)r+[^ + ^)]|5t(q)| 2 

IML, q L L Po ^lPo Po 

--[^(q)M-q) + ^(-q)Mq)]) + ^ > (s) 

Po J 

where F' is the equilibrium value. The Gaussian approximation to our hy- 
drodynamic model for the pure system can be viewed as a modified Ornstein- 
Zernike theory for flux lines [29]. The familiar Gaussian free energy of di- 
rected lines is obtained by setting Sp(q) = St z (q) in Eq. (5) and neglecting 
terms of higher order in q, with the result 

^ [6t] * F ' + ^LpJ ^ {^(^Uq)! 2 + c?i(?)l^(q)l 2 } + ?d , (6) 

where 

c° u (q) = Po ei + p 2 V(q) , (7) 
<%M=f%V(q), (8) 

are the bare local tilt and compressional moduli, respectively, of a pure di- 
rected flux-line liquid. 

The correlation functions of the fluctuations of the hydrodynamic fields 
are immediately evaluated by carrying out both a thermal average with 
weight ~ exp[— jF G //c B T] and an average over quenched disorder. Retain- 
ing here and below only the first-order corrections in A, we obtain 

(9) 



PoQl 



(^(q)^(-q)) G = ^T^ + A 



PoQl 
D(a). 
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(^(q)M-q)>o = k B r 



+A 



£>(q)" 



^(q) 



^(q±) , 



(10) 



(5p(q)5p(-q)) G = A^T 



+A 



Po<? 2 c^ 4 (g) 



£>(q) p ei(l + g 2 ag(a /2/ 2 ) 2 ) 



Po9. 



2 



£>(q)(l + g 2 a 2 l (a /2Z 2 ) 2 ) 



2 



(11) 



where i, j = (x,y), i^(q±) = q±iq±j/ql and i^(q±) = % - ^(q±) are 
longitudinal and transverse projection operators and 

£(q) = q 2 AM/Po + ql[c°n(q) + <W?)]/>o • (12) 

As seen below, <5cn(g) represents a short scale correction to the compressional 
modulus. It is given by 

r / \ q 2 a 2 (a /2l z ) 2 

5Cll(9)=Poei l + «(a /2/ z ) 2 ' (13) 

with l z = eial/(2k B T). The autocorrelation function of the longitudinal part 
of tj_ is simply related to that of t z as the "continuity" equation (2) requires 
q_i_ • tj_ = —q z t z . Cross correlations among the fields are also not vanishing 
but will not be given here. 

Compressional and tilt deformations of the line liquid can be described 
by defining effective local (wavector-dependent) elastic moduli as 

p 2 k B T 



c u (q±) Qz-*0 
p\k B T 



lim(<ft 2 (q)<M-q)>, (14) 



cu(fl.) = ^^ (qi)<t|(q) ' i( - q)> - (15) 

In the absence of disorder, this gives Cn(q±) = c? 1 (gj_) + 5cu(q±), while 
cu(q±) = c® 4 (q±) is unchanged. The length scale l z is the thermally-induced 
entanglement correlation length introduced by Nelson [12]. A single flux line 
at finite temperature wanders transversally as it crosses the sample along 
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the z direction. The entanglement length l z represents a "collision" length 
in the time-like variable z. It is the length scale a single flux line travels 
in the z direction before wandering a transverse distance of order a and 
colliding with a neighbor. Such "collisions" tend to stiffen the compressional 
modulus on small length scales (of order a ) due to the caging of each line 
by its neighbors. 

The long-wavelength disorder-renormalized elastic constants are defined 

as 

p 2 k B T 



cfi '/ -(>,/-.-() 



lim .Urn (8t z (ci)8t z (-q)) , (16) 



p 2 k B T 



^f- = lim 1m p£(cu)<Mq)*,-(-q)> , (17) 

where the order of the limits is important. These definitions are consistent 
with the definition of the elastic constants as response functions. Within the 
Gaussian approximation, point disorder described as a random field coupled 
to t z (r) softens the compressional modulus, but it does not renormalize the 
tilt modulus, 



^ l + Apl/(cUO)k B T) 



(18) 



C44 = 4(0) • (19) 

As discussed in Ref. [26], the long- wavelength limit of the autocorrelation 
function of the fluctuations in the field t z (r), associated with spontaneous 
vortex loops, determines the renormalized single-line energy ei, according to 

hm lim MWFl^ = ^ , (20) 

irrespective of the order of limits. A random field that couples to t z gives 
no corrections to the line energy and ef = e\. It is important to stress the 
difference between the line energy, ei, and the line tilt stiffness, £\. The former 
is the energy per unit length needed to create a flux line in the direction of 
the external applied field. It is related to the lower critical field and, in the 
boson analogy, it corresponds to the boson rest mass. The tilt coefficient 
e~i measures the stiffness of a line against a tilt away from the direction of 
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the external field. In clean, isotropic superconductors, e\ = e~i are equal, but 
they can differ strongly in anisotropic materials [30]. Using the free energy of 
strongly fluctuating flux-line liquids in uniaxial anisotropic superconductors 
[26] (where the external field is applied along the c-axis of the material), one 
can show that the long-wavelength limit of autocorrelator of the t% field yields 
ei, while the single- vortex part of the tilt modulus yields e~\. In the following 
sections, we shall show that different types of disorder have different effect 
on these two distinct quantities. 

Finally, the disorder-induced correction to the autocorrelators of the two 
different densities, Eqs. (9) and (11), differ at finite wavevectors. This dif- 
ference becomes pronounced in strongly entangled liquids (high temperature 
and/or density), where ao/l z ~ 1 (a /l z ps 0.3 for T = 77K and B = 2T). 

4 Scalar point defects 

In this section we model uncorrelated point disorder as a scalar random field 
that couples to the scalar field p(r), 



This corresponds to a random contribution to the line stiffness e\ (or to 
the quantum particles' mass in the relativistic charged boson mapping). The 
random field is assumed to have zero mean and correlations given by Eq. (??), 
although with a different disorder strength, A. This random field physically 
represents uncorrelated pinning centers which favor the nucleation of flux 
lines irrespective of their orientation. We therefore expect that it will favor 
vortex loop formation. In a liquid of directed lines, there is no difference 
between this case and the one discussed in Section 3. 

The thermal part of all correlation functions is, of course, unchanged. 
Introducing the notation S{...) = (...) — (...) to denote the contribution due 
to disorder to the various correlation functions, we find 




(21) 



5(5t z (q)5t z (-q)) G 



A 



(22) 



D(q)(l + q 2 al(a /2l z y) 




(23) 
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5(5p(<i)5p(-q)) G = A 



Po<f c° M (q) 



(24) 



_D(q) p ei(l + q 2 al{a /2l z 

Again, disorder leaves the transverse part of the tj_ correlations un- 
changed and therefore it does not renormalize the long-wavelength tilt mod- 
ulus (in the Gaussian approximation). The renormalized long-wavelength 
compressional modulus is identical to that given in Eq. (18), although the 
two differ at finite wavelengths. This type of scalar disorder does, however, 
yield a finite renormalization of the line energy that determines vortex loop 
fluctuations. Using the definition given in Eq.(20), we find 

Ap 

The disorder-induced correction to the line energy on the right hand side of 
Eq. (25) can be written as the ratio of two length scales as ef = ei/(l + 
Iz/Id), where l z is the thermal entanglement length introduced earlier and 
Id = e^/(2Apg) is a characteristic entanglement length associated with dis- 
order. This type of quenched disorder softens the line tension, enhancing 
the formation of vortex loops. The condition l z l D marks a crossover at 
T ~ (A/ei)B z from a regime where loop fluctuations are thermally domi- 
nated to one where they are induced by quenched-disorder. Although nu- 
merical studies of the pure uniformly frustrated XY model have ruled out 
the scenario of the "loop blow-out" transition proposed by Tesanovic [20, 21] 
and by Sudb0 and collaborators [22], our result indicates that vortex loop 
formation is enhanced by disorder, suggesting the possibility that such a 
transition may be recovered in sufficiently disordered systems, where the line 
tension may be driven to zero. 



5 Point defects with randomly oriented mag- 
netic moments 

Finally, in this section we consider uncorrelated point impurities that can 
be modeled as quenched magnetic dipoles with random orientation. This 
corresponds to a vector disorder potential that couples with statistically equal 
strength to all components of the local magnetic field, yielding 

Td = Jd 2 r ± dzV D (r)-t(r), (26) 
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with V DlI (r) = and V Dii (y)V Dv (y') = 5^ u A5(r - r'), for fi,u = (x, y, z). The 
random field favors alignment of the local magnetic field due to flux lines 
with uncorrelated quenched random directions. 

To Gaussian order, the disorder-induced corrections to the correlation 
functions are 



5(5t z (q)5t z (-q))G = A 



(27) 



*<fi(q)fi(-q)>G = A 



P% 



^(q±) + 



2 ^(q±)}, (28) 



8(8p(q)8p(-q)) G = A 



poQ±q 



D(q)(l + q*al(a /2l z )i) 



(29) 



The disorder-induced correction to the oriented areal density autocorrela- 
tion, Eq (27), differs significantly from that obtained for the other two types 
of disorder when q± pa and q z ^ 0. 

As in the case presented in Section 3, this type of disorder does not 
renormalize the line stiffness e\. The transverse part of the tangent field 
autocorrelator does, however, acquire a finite correction, which leads to a 
renormalization of the long- wavelength tilt modulus. From Eq. (17), we find 



1 1 A 1 

~ & + k B T [ ei + p o 0g/(47r)f 



(30) 



As shown by Larkin and Vinokur [31], the renormalization of C44 can be 
related to that of its single- vortex part, c v 4A = po^i, with 1\ the tilt stiffness 
of a directed line. This is renormalized by the "vector" disorder according to 



ef = e 1 (l - l z /l D ) . 



(31) 



Isotropic random dipoles tend to align flux-line segments along the random 
directions of the magnetic moments carried by the impurities and facilitate 
tilting away from the z-direction. On the other hand, the line tension defined 
in Eq. (20) and associated with loop fluctuations remains unrenormalized 
by this type of disorder. In both cases, the disorder- induced corrections 
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are determined by the ratio I z /Id of the thermal and disorder entanglement 
lengths (see Eq. (25)). 

Finally, the disorder-induced renormalization of the compressional mod- 
ulus is the same as for the other two types of point disorder in the long- 
wavelength limit, although it differs at finite wavevectors. 

6 Conclusions 

In this paper, we have calculated weak- disorder- induced corrections to the 
static correlation functions of strongly fluctuating (arbitrarily curved and 
looping) flux-line liquids using a Gaussian hydrodynamic approach. Three 
realizations of point disorder have been considered, corresponding to mag- 
netic and non-magnetic defects. Our results for the structure function, the 
scalar density (p(r)) autocorrelator, and the tilt field (tj_(r)) correlations 
can, in principle, be probed experimentally using polarized and unpolarized 
neutron scattering. We have predicted qualitative differences between these 
quantities at finite wavevectors for different types of disorder which are ex- 
perimentally testable. From these correlations, we have also obtained the 
disorder-induced corrections to the elastic moduli. The compressional mod- 
ulus gets softened for all three types of point disorder. Loop fluctuations 
are enhanced by scalar disorder. The tilt modulus is softened by isotropic 
random polar (magnetic) disorder. 

A softening of the elastic constants by quenched disorder may appear 
surprising, as naively one would expect pinning to stiffen the system. A 
result qualitatively similar to ours has, however, been obtained before in a 
different system [32]. Preliminary results of a calculation of the disorder- 
induced correction to the compressional modulus taking into account non- 
Gaussian terms in the free energy suggest that this correction may change 
sign depending on the system parameters (temperature or magnetic field), 
indicating that a perturbative approach may not adequate in this case [33]. If 
the softening were to survive a non-perturbative calculation, this would signal 
a disorder-driven instability of the vortex liquid phase. Clearly more work 
beyond the Gaussian approximation is needed to answer these questions. 

MCM was supported by the National Science Foundation under grant 
DMR-0305407. 
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